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. . . . ABSTRACT 
O 

o 

^ ' We calculate the coupled thermal evolution and magnetic field decay in rela- 

■ tivistic model neutron stars threaded by superstrong magnetic fields {B > 10^^ 
^ ■ G). Our main goal is to evaluate how such "magnetars" evolve with time and how 

■ field decay modifies the transitions to core superfiuidity and cooling dominated 
^ ■ by surface X-ray emission. Observations of a thermal X-ray spectral component 

and fast timing noise place strong constraints on the presence of a superfiuid 
core. We find that the transition to core superfiuidity can be significantly de- 
^ ■ layed by field decay in the age range ~ 10^ — 10^ yrs. The mechanism of Hall 

drift is related to the stability of the core magnetic field, and to currents fiowing 
^ ' outward through the crust. The heating effect is enhanced if it is continuous 

O ■ rather than spasmodic. Condensation of a heavy element layer at the surface is 

• shown to cause only modest changes in the outward conduction of heat. 

2 ' Subject headings: magnetic fields: stars - neutron stars: general 

> 

X 

' 1. Introduction 

Observations of the Soft Gamma Repeaters and Anomalous X-ray Pulsars have revealed 
a rich phenomenology of transient X-ray emission and torque variability. Some members of 
the first group occasionally emit enormously bright X-ray fiares (Hurley 2000), which are 
followed by transient periods of decaying X-ray fiux. The unified nature of the SGRs and 
AXPs is indicated by the detection of ~ 0.1 s hard-spectrum X-ray bursts from two AXPs 
(Gavriil et al. 2002; Kaspi et al. 2003). This variable X-ray emission is now generally 
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believed to be powered by the decay of an ultrastrong magnetic field (Thompson & Duncan 
1996; Colpi, Geppert & Page 2000). 

Two sub-classes of these objects have emerged recently: those which maintain a fairly 
steady X-ray luminosity over two decades or longer; and a less conspicuous group (not 
necessarily smaller in total numbers) in which one sees transitions to and from very low 
levels of persistent X-ray emission over a period of months to years (Torii et al. 1998; 
Kouveliotou et al. 2003; Ibrahim et al. 2003). Objects in the second group are weak burst 
sources or have not been observed to burst at all. 

Some SGRs and AXPs also show strong timing irregularities, which broadly can be 
divided into i) adiabatic variations in torque over a period of months to years (Kaspi et al. 
2001; Woods et al. 2002) and ii) a broad spectrum of timing noise starting from periods as 
short as ~ 10"^ s and upward. Our goal here is to isolate the most important mechanisms by 
which an ultrastrong magnetic field will contribute to the persistent X-ray emission of these 
objects (on timescales of years or longer); and also to make deductions about the composition 
of the star from the appearance of strong timing noise in active burst sources. 



2. Basic Modes of Magnetic Field Decay 

A decaying magnetic field in a neutron star evolves through a series of equilibrium states, 
punctuated by the release of elastic stresses in its crust and the excitation of hydrodynamic 
motions in its liquid core. The crust is too weak to sustain all but modest departures from 
magnetostatic equihbrium, ii B > 10^^ G. In the case of a purely fluid star that is stabihzed 
against convection in the radial direction, it is generally believed that the equilibrium states 
of the magnetic field must carry some net helicity, e.g., that the field has both toroidal and 
poloidal components. Purely poloidal (Flowers and Ruderman 1975) and purely toroidal 
(Tayler 1973) fields are unstable. 

The magnetic field evolves according to the equation 

= V X [(V + Vamb + VHall) X B] + — (l) 
Ut Ot fracture 

Here Vamb is the speed with which the magnetic field is advected by the diffusing charged 
component (electrons/protons) of the neutron star core; VHaii = — J/erZe; and v is the hy- 
drodynamic response of the core to the combined effect of these transport processes. The 
rigid crust is also subject to sporadic yields and fractures which cause changes in B (both in 
crust and core) on short timescales. The star is spherical in a first approximation, and the 
compositional stratification enforces Vj. — 0{eB) — ^ and V • v = 0(£b) ~ (Reisenegger & 
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Goldreich 1992). Here Eb — B'^/SnPe and Pg is the electron pressure. 

The diffusion of charged particles in a normal, degenerate n-p-e plasma is limited by 
proton-neutron drag at high temperatures (Haensel et al. 1990; Goldreich & Reisenegger 
1992, hereafter GR); and at T < 5 x 10^ K by the rate of relaxation of chemical potential 
gradients (GR; Pethick 1992). We note here that the n-p coUision rate is also reduced 
dramatically after the transition to proton superconductivity (due to the lower density of 
quasiparticle excitations near the proton fermi surface). This probably happens early on 
compared with the ~ 10'' yr lifetime of SGR/AXP activity (e.g. Yakovlev et al. 2001). To 
calculate the continuing ambipolar diffusion of a magnetic field in a superconducting core, 
we therefore include only the limiting effects of chemical potential gradients. In this case, the 
timescale is tamb(-B,T) ~ SirnKkBTy /UTj^ca.{p,T)B'^ when the chemical potential imbalance 
AfjL ~ B'^/SnUe induced by the J x B/c force is much less than ksT. Here C/urca is the rate 
at which beta transformations between neutrons and protons release energy to neutrinos. 
This rate is strongly modified by Cooper pairing of protons and neutrons, and we include the 
resulting corrections as tabulated in Yakovlev et al. (2001). For mean fields stronger than 
several xlO^^ G, the magnetic sheaths of the superconducting fiuxoids are packed together. 
Any collective motion of the fiuxoids with respect to the protons is strongly inhibited, and 
we can treat the field as being continuous. In a normal n-p-e core, the heating induced by 
the decay of the field leads to a direction relation ksT ~ A/x oc B'^ between T and B. Given 
the strong T-dependence of the URCA rates, this causes a strong feedback iamb on 
the drift rate above a critical flux density of ~ 3 x 10^^ G (Thompson & Duncan 1996). 

If the interior magnetic field of a neutron star is helical, then the current density J and 
B both have poloidal components. In a magnetar, closure of this poloidal current inside the 
star generates a J x B/c force which is strong enough to fracture the crust, thereby twisting 
up the external magnetic field (Thompson, Lyutikov, & Kulkarni 2002). This mechanism 
is driven by Hall effect in the crust, through the gradient in electron density with height 
z. The Hall term in eq. (1) yields t^^^ = B~^^{dB^/ dt) = Jz{d / dz){ene)^^ in cylindrical 
coordinates. Where the poloidal current flows outward - > - the toroidal field will 
migrate toward lower densities. There the crust has a smaller shear modulus and is less able 
to balance magnetic shear stresses. The net effect is to excavate the twist in the poloidal field 
from the crust, on the timescale tnaii = 2.4 x 10^ {Rq/ B^^i^)p^^l yr at a radius R — RqX 10^ 
cm and density x 10^^ g cm"^. 

The Hall effect has different consequences in the fluid core. In a liquid that is stratifled 
parallel to gravity g, Hall drift creates unbalanced magnetic stresses. In particular, hall 
waves do not exist as propagating modes in the high frequency limit, k\d\n p / dz\^^ ^ 1. To 
show this, note that the Hall term can only be cancelled by a hydrodynamic displacement 
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satisfying Vz — 0. This condition is satisfied if Vnaii has no vertical component, i.e., if = 0. 
A simple example is a finite-amplitude field variation B' = BQexp(ik • x) superimposed 
on a uniform background field Bq. Magnetostatic equilibrium requires that B' be linearly 
polarized, and that (Bq x k) • g = 0, which implies in turn — 0. 

On larger scales, the Hall term in the induction equation generally cannothe cancelled off 
by such a hydrodynamic displacement field in the fluid corc.^ Hall drift evolves the magnetic 
field into a new configuration of equal total energy (GR), but if the field is initially an a 
stable configuration then all displacements satisfying = 0, V ■ v = increase its energy. 
In fact, the new configuration will generally have a higher energy than some neighboring 
configuration into which it can relax. In what follows, we make the reasonable assumption 
that this excess energy is converted to heat on a timescale much less than ~ 10^ yrs. 

The next step is a prescription for the hydrodynamic response of the fiuid interior to 
this diffusive motion. We define a characteristic poloidal magnetic field Bp, toroidal field 
Bij) = O{10)Bp, and tilt angle a between the poloidal current and the direction of gravity 
at the center of each toroidal loop. There is a characteristic tilt, a ~ Bp/ B^, below which 
the field is able to largely unwind through differential rotations of fluid shells conflned to 
gravitational equipotential surfaces (Fig. 2 of Thompson & Duncan 2001). Hence changes 
in a can cause unwinding of the toroidal field at the rate dB^p/dt — —(B^/a) \da/dt\. On 
the other hand, microscopic transport of the field on a timescale r causes changes in a at a 
rate da/dt ±r~^. Thus we evolve the simple one-field model 



Unwinding of the field does not require any change in Bp, and the associated transport time 
is much longer when Bp <^ B^. Hence we only evolve B^. 

A net change in winding of the core magnetic field must be accompanied by a torsional 
deformation of the stellar crust. The energy deposited in an area A of the crust by twisting 
it through an angle 9 is 9^ J fxdV ~ 5 x 10^° {9/0.001)^ (A/100 km^) ergs. (Here = 1.1 x 
10^°p°f is the crustal shear modulus; Strohmayer et al. 1991.) During a concentrated episode 
of SGR activity, when hundreds of ~ 0.1-s X-ray bursts are emitted, this process might be 
repeated ~ times, and the net crustal heat deposition would be ~ 1000 (^/O.OOl)"^ times 
larger. A giant flare evidently involves a single readjustment through a much larger angle. 



-'^Mestel (1956) has argued that Hall drift of an axisymmetric poloidal magnetic field will excite a torsional 
Alfven wave in a fluid star. However, such a field configuration has a special symmetry (radial current 
Jr = 0) and the above argument indicates that the Hall drift will instead be cancelled by a compensating 
hydrodynamic motion. 



dB^ 



B^ \ 1 ^ 1 



(2) 



dt 



Q; L^amb(-B</,) tHaXliBtf,) 
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A toroidal magnetic field confined to the crust will, in a first approximation, evolve 
independently from that in the core. The net angle through which a radial poloidal field 
is twisted (across one vertical density scale height) is limited by the finite yield strain i/j 
of the crust, A0b = {ip/ R){B^/ B^) < {Anfi/ B^)iIj. This sets an upper bound B^ < 7 x 
10^^ B~1^{iP/10~^)Rq G. Integrating the above Hall equation for B^p, under the assumption 
that inhomogeneities in the winding of the field relax completely in the deep crust after 
each yielding event, gives B^/ B^^q = [1 + t/tHaii(-B</,,o)]~^- The corresponding heating rate, 
integrated over depth 2 x 10^^ j^gAiai^g ^l+i/iHaii(-B0,o)]~^ ergs s~^, is constant as i — > 
and scales as at i ^ ^Haii(-B</,,o)- We will revisit the interplay between radial currents and 
Hall drift elsewhere in more detail. 



3. Coupled Thermal Evolution and Field Decay 

The thermal evolution of the star is followed using a relativistic model constructed 
from the Tolman-Oppenheimer-Volkoff equations and the parametrized equation of state of 
Prakash, Lattimer, and Ainsworth (1988) with intermediate compressibility {K — 240 MeV) 
and a mass low enough that direct URCA cooling can be neglected (M = 1.35 Mq). (We note 
that more recent EOS indicate the absense of direct URCA cooling for much larger masses; 
Akmal, Pandharipande, and Ravenhall 1998.) We are interested here in baseline thermal 
evolution, averaged over a timescale long compared with the conduction times across the 
crust, which is more than two orders of magnitude shorter than the spindown ages of the 
SGRs and AXPs (P/2P > 10^ yr). The calculation therefore takes into account the effects 
of continuous modes of magnetic field decay. 

We calculate the net neutrino luminosity by integrating the emissivities due to the 
modified URCA reaction ({n,p} + n ^ {n,p}+p + e~ + i/e), nucleon bremsstrahling emission 
(n+n — > n+n+p+i', etc.) and neutron Cooper pair emission {n-\-n —>■ (2n)+i/+i?); evolving 
the magnetic field according to eq. (2) using the volume- averaged modified-URCA rate in 
^amb; calculating the luminosity Lx — 'inRljge'^'^'^^^^^asBTg (0(r) = gravitational potential) 
in thermal surface X-ray emission from the relation between core temperature and surface 
temperature Tg detailed below; and, finally, evolving the (rcdshifted) core temperature Tc = 
e'^^^^T{r) in the isothermal approximation by balancing the net cooling luminosity with the 
rate of loss of magnetic energy from the volume V of the star. 

Balancing only the last term with the modified-URCA emissivity, one finds A// ~ B'^/SnUe ~ 
a{kBTc). There are finite- A// corrections to the URCA rates (Reisenegger 1995), but these 
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are generally less important than the corrections due to nucleon pairing when a <^ 1. 

The heat flux emerging through the surface of a magnetar depends on the relation 
between the surface effective temperature Ts and the temperature Tc in the deep crust. 
This relation has been calculated (Van Riper 1988; Heyl & Hernquist 1997; Potekhin and 
Yakovlev 2001) for a magnetized atmosphere in which the ions form a non-ideal gas. At the 
surface temperature characteristic of SGRs and AXPs (fc^Tbb ~ 0.4-0.5 keV; Ozel, Psaltis, & 
Kaspi 2001) one however expects heavy ions (e.g. iron) to be condensed into long molecular 
chains (Lai and Salpeter 1997). The density is then large even close to the surface and can 
be estimated by minimizing the sum of the Coulomb and electron degeneracy energies in a 
Wigner-Scitz ceU:^ p(0) ~ 1.8 x W {YjQ.b)-\Z /2&f/^{B /W^ Gf/^ g cm'^. 

The transmission of heat through the atmosphere of a neutron star with zero surface 
magnetic field is controlled by a sensitivity strip of a lower density p ~ 10^ — 10^ g cm~^ 
(e.g. Potekhin & Yakovlev 2001). Fig. 1 shows the surface temperature of a neutron star 
with normal core neutrons and protons (transition temperatures Tc,„ = Tc^p — 0) using 
this standard envelope model. We have also included the effect of ambipolar diffusion of 
an internal toroidal field (S^ = 3, 5 x 10^^ G), using eq. (2) and the ambipolar diffusion 
timescale iamb tabulated in eqs. (58), (59) of GR. (See also Heyl & Kulkarni 1998.) 

The heat flux is higher in a second envelope model which we have constructed for a 
condensed iron surface layer in a 10^^ G magnetic field. The best power-law fit is Tg/ g\^^ — 
4.0 X 10^ K {Tjlif K) V2 [2.9 x 10^ K {Tjlif K) V^] for greater than [less than] 2.8 x lO^K, 
and gravity g — gu x 10^^ cm s~^. The integration of dT/dz over depth z ideahzes the 
electrons as a (locally) uniform fermi gas, and employs the thermal conductivity tabulated by 
Potekhin (1999). The ideal finite-T electron equation of state is supplemented by a correction 
for the (negative) Coulomb pressure in the solid phase, P = P^. — Pf,{0)[ne/ne{0)]^^^ , where 
Pe(0) and ne(0) are the electron degeneracy pressure and density at zero total pressure P. 

Much of the thermal resistance is localized where the electrons are 1-dimensional, with 
fermi momentum ppe ~ "^eC (density p = YempeBppe/^Ti'^li^c ~ lO^-Bis g cm^^) and where 
the melting temperature is 1 x 10^(^/26)^/^5^5^ K. Thus at high the sensitivity strip lies 
in a Coulomb liquid layer below the solid surface; whereas at low Tc (and large Z) it lies 
within the sohd. In the sohd, the electrons are nearly degenerate and heat is transported 
by conduction electrons within a narrow energy range ~ ksT near the fermi surface. The 
scalings for Ts(Tc) are easily derived by approximating the conductivity as being due to 
phonon scattering (in the solid) or Coulomb scattering (in the liquid) , with all the electrons 
in the lowest Landau state. This analytic model gives oc bI^^Z~^/^ at high Tc, and 



^Here Ye is the electron fraction and Z the nuclear charge. 
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Tg oc Bg at low Tc. To correct for the lower opacity of the extraordinary mode near the 
condensed surface, we show 2^/^ times the effective temperature. Since there is evidence for 
complicated multipolar structure in the surface magnetic field of some SGRs (Feroci et al. 
2001), Bs is a characteristic surface field without any dipole structure. 

Our numerical results are plotted in Figs. 1 and 2. In sum: 

i) Steady heating of the stellar interior by a ~ 3 x 10^^ G magnetic field increases 
logio(T's) by ~ 0.2 if the core neutrons and protons are normal. Proton pairing suppresses 
the rate of charged-current weak interactions and therefore the rate of the irrotational mode 
of ambipolar diffusion; in which case Hall drift makes a similar contribution to the decay of 
the core magnetic field before the core superfluid transition, and dominates thereafter. 

ii) If the peak pairing temperature Tc^„ ^ 6 x 10^ K, the drop in Tg is gradual and 
occurs earher than ~ 100 years; but if Tc,„ < 5 x 10^ K, then this drop is sharper and occurs 
in the observed age range of SGR/AXP activity. Magnetic dissipation can delay the time 
of the pairing transition by an order of magnitude. It also delays (and makes sharper) the 
transition to photon-dominated cooling. If the internal B-ficld is much stronger than 10^"^ 
G, and its coherence length is large (~ 10 km), then continuous Hall decay (which is not 
temperature sensitive) will keep the surface warm up to an age of several xlO^ yrs (redshifted 
temperature 7;°° = T^e-^^^Ns) 1.5 x 10*^ K). Strong intermittency in the rate of Hall decay 
would allow the core to cool off more rapidly in between magnetic 'flares', due to the strong 
temperature dependence of the neutrino emissitivities, and the temperature evolution would 
be intermediate between the dashed and sohd curves in Fig. 2. 

iii) The thermal evolution with a condensed Fe atmosphere gives 7"^°° within a factor 
~ 1.3-1.5 of those measured in the thermal components of AXP spectra (/cTbb — 0.4-0.5 keV) 
at an age of ~ 10^-10^ yrs. A relatively thin light-element layer will harden the thermal peak 
by this amount, compared with a pure black body (e.g. Lloyd, Hernquist, and Heyl 2003). 
If the surface field is ~ 10^^ G, then Lx is generally less than 10^^ erg s~^, consistent with 
some AXP sources but not all (e.g. Ozel et al. 2001). Recalculating the cooling models 
with with a lighter element (carbon) surface layer shows a factor ~ 3 — 4 increase in thermal 
transparency at early times - but not below ~ 4 x 10® K, so that the photon cooling time 
does not change significantly. Heating by an external current may also not be negligible: the 
power needed to force a current of ions and electrons through a twisted magnetosphere with 
a constant equatorial pitch B^/Bq is Lx ~ 3 x 10^^ (Spde/lO^^ G){B^/Bq) ergs/s (Thompson 
et al. 2002). 
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4. Timing Noise and Core Neutron Superfluidity 

Phase-resolved X-ray timing of SGRs 1806—20 and 1900+14 have revealed a broad spec- 
trum of timing noise (Woods et al. 2002). At short time invervals {At ~ 3 x 10^ s) large 
phase offsets (A0 ~ 0.3 cycles) are sometimes observed in the X-ray pulses, which prevent 
phase-connected timing. These offsets correspond to stochastic shifts in frequency of mag- 
nitude Aiy/u ~ 10~^. The alternative is an enormous increase in frequency derivative by 
~ 10^ (Az> ~ 2A0/Ai^ ~ 10~^ s~^). Hence the effect appears to be due to the absorption 
and release of angular momentum by an internal superfluid component. The lack of obvious 
classical post-glitch relaxation behavior, and the magnitude of the timing residuals, lead us 
identify a superfluid core as the reservoir of angular momentum. Indeed, no obvious corre- 
lation exists between this timing noise and the observation of long-term torque variations, 
bright X-ray outbursts, or large changes in persistent X-ray flux that would be associated 
with deformations of the crust. 

Timing noise of this amplitude is absent in radio pulsars, and its presence in magnetars 
is strongly correlated with overall activity as a burst source, and with a spectrally hard 
persistent X-ray emission. Part of this correlation could result from the reduction in interior 
temperature (by a factor ~ 5) coinciding with the pairing transition of the neutrons, which 
would make the crust more brittle; and from a reduction in surface cooling rate in comparison 
with the power dissipated by external currents. 

In a slowly rotating magnetar, the outward motion of the vortices is constrained by the 
magnetic field, due to the large energetic barrier to the crossing of vortices and supercon- 
ducting fiuxoids (e.g. Ruderman et al. 1998). Slow fluctuations in the field of amphtude 
AB/B ~ [Isi/ I)'~^Av /v ~ 10~^(/sf/-?^)~^, with some amplitude perpendicular to the axis of 
rotation, would provide the necessary perturbation to the superfluid. If such fluctuations 
are to occur over many (> 10*^) core Alfven-crossing times, the field configuration must be 
metastable, and nearly degenerate in energy with others that differ by one part in ~ 10^. 
This is expected if the core field has relaxed significantly from an initial equilibrium state. 

5. Conclusions and Observational Tests 

Several degrees of freedom strongly influence the thermal and spectral evolution of 
a magnetar: i) the strength of its internal magnetic field; ii) the temperature at which 
core neutrons become superfiuid; iii) the transmissivity of its thermal envelope; and iv) the 
configuration of its internal magnetic field and the equilibrium states through which this 
field moves. Before a core superfiuid transition, thermal luminosities up to ~ 10^^ ergs/s 
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are possible if the surface magnetic field is ~ 10^^ G; thereafter, the neutrino emissivity 
rises dramatically due to Cooper pair cooling, and thermal emission powered by internal 
Hall decay can continue at a level < 10'^^ ergs/s beyond an age of ~ 10^ yrs. Wc infer 
that magnetar candidates which show transitions to/from such low luminosity states have 
magnetic fields of intermediate strength and/or superfiuid cores. We have also argued that 
the observation of fast timing noise in SGRs 1900+14 and 1806—20 provides direct evidence 
for a superfiuid neutron core. The highly non-thermal persistent X-ray emission of these 
sources must be powered mainly by external currents; but a thermal seed could be provided 
by internal heating. (The AXP with the hardest X-ray spectrum and noisiest spindown, IE 
1048.1—5937, also has a relatively low blackbody luminosity, less than 10^^ ergs/s; Ozel et 
al. 2003.) SGR 0526—66 has remained X-ray bright since its last observed outburst in 1983 
(Kulkarni et al. 2003) and has a high luminosity ~ 7 x 10'^^ erg/s: either B ^ 10^^ G at its 
surface, or most of its emission continues to be powered by external currents (in spite of the 
relatively soft X-ray spectrum). Detailed timing measurements of this source would provide 
valuable diagnostics of its interior state. 

We thank Lars Bildsten, Vicky Kaspi, Shri Kulkarni, Bennett Link, and Peter Woods 
for conversations. PA is an AAPF NSF fellow, AG is a Hubble Fellow, and GT acknowledges 
the support of the NSERG of Ganada. 
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Fig. 1. — Surface temperature (redshifted) as a function of age, for a 1.35Mq neutron star 
with normal neutrons and protons. Upper two curves in each set correspond to stars with 
internal heating mainly by ambipolar diffusion.) Bg is the surface field, which is assumed 
to be constant in time, and the internal toroidal field (which may be stronger). See text 
for the description of the envelope models. The factor of 2^/^ corrects for surface emission 
dominated by the extraordinary polarization mode. 
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Fig. 2. — Same but for a condensed iron envelope, superconducting core protons (Tc,p = 
5 X 10^ K), and various peak pairing temperatures for the core neutrons. The distribution of 
critical temperature Tc^„ with density is taken to be Tc^„ = Tc^„ (peak) exp[—(pi4 — 5)^/p^4]. 
Smooth declines in X-ray flux at t > 10^~^ yrs occur when surface cooling begins to dominate. 
Sharper drops signal a transition to core superfluidity. Those occuring at i < 100 yrs 
correspond to rc,„(peak) = 9, 7 x 10^ K; and at t ~ 10^ yrs to rc,„(peak) = 5 x 10^ K. Sohd 
curves {Bg — 10^^ G) include the effects of Hall decay. Dashed curves do not, and ordered 
from left to right correspond to rc,n(peak) = 9, 7, 5, x 10^ K. Notice that the superfluid 
transition is greatly delayed for Tc^„ = 5 x 10^ K. Dotted curves show weaker magnetic fields, 
= fi^(O) = 2 - 10 X 10^^ G and Tc,„ = 5 x 10^ K. Finally, the top dashed curve shows 
enhanced cooling through a condensed carbon envelope (density < 1 x 10^° g cm~^) with 
B, = lx 10^^ G, B^{Q) = 3 X 10^^ G, and Tc,„ = 5 x 10^ K. 



